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Abstract 

In this paper we discuss the theoretical difficulties in extracting Vub using the 
data from inclusive B decays. Specifically, we address the issue of the end point 
singularities. We perform the resummation of both the leading and next to 
leading end point logs and include the leading corrections to the hard scattering 
amplitude. We find that the resummation is a 20% — 50% effect in the end point 
region where the resummation is valid. Furthermore, the resummed sub-leading 
logs dominate the resummed double logs. The consequences of this result for a 
model independent extraction of the mixing angle Vub are explored. 



1 Introduction 



Measurements in the bottom quark sector have reached the point that our knowledge of many 
observables is now bounded by the theoretical uncertainties Fortunately, theoretical 
advances in calculating both exclusive as well as inclusive rates now allow the extraction of the 
CKM parameters without recourse to the models which have soiled the extraction processes 
to data. The present values of Vub have a model dependence which introduce an uncertainty 
of a factor of 2[Q, which is several times larger than the experimental uncertainties. With 
QCD based calculations, we can now hope to extract both Vbc and Vub with errors on the 
order of tens of percents. In this work, we concentrate on the extraction of Vub from the 
measurement of the electron spectrum in semi-leptonic inclusive B meson decays. 

The extraction of Vub from inclusive semi-leptonic B decays is hindered by the fact that 
the background from charmed decays is overwhelming for most of the range of the lepton 
energy. Thus, we are forced to make a cut on the lepton energy, vetoing all events, or some 
large fraction thereof, with lepton energy less than the 6 — >■ c end point energy. Given the 
proximity of the two relevant end points, this obviously hinders the statistics. However, 
even with a large data sample, the accuracy of the extraction will be limited by the errors 
induced from the approximations used in calculating the theoretical prediction in the end 
point region. This region of the Dalitz plot is especially nettlesome for theory, because 
the perturbative, as well as the non-perturbative corrections become large when the lepton 
energy is near its endpoint value. 

It has been shown that it is possible to calculate the decay spectrum of inclusive heavy 
meson decay in a systematic expansion in e = and as using an operator product expan- 
sion within the confines of heavy quark effective field theory^. It is possible to Euclideanize 
the calculation of the rate for most of the region of the Dalitz plot with only minimal as- 
sumptions about local duality. However, in the end point region, the expansion in e, as well 
as the expansion in a.^, begin to breakdown (The endpoint region poses problems for local 
duality as well. We shall discuss this in more detail later). 

The aim of this paper is to determine the size of the errors induced from the theoretical 
uncertainties in extraction of Vub- A large piece of this work consists of implementing the 
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resummation of the leading and sub-leading endpoint logs which cause the breakdown of 
the expansion in a^, as first discussed on general grounds in ||^, and the inclusion of the 
Us corrections to the hard scattering amplitude. However, to determine the consistency of 
our calculation, we must also address the issue of the non-perturbative corrections. These 
issues have been previously looked at in refs 0] and 0. In [Q the need for resummation 
was addressed on general grounds. However, the calculational methods used here are not 
compatible with the arguments given in 0], and thus we must recapitulate these arguments 
within the confines of our methods. 

In the second section of this paper, we discuss the question of the need to resum the 
perturbative as well as non-perturbative series. The next three sections are dedicated to 
the resummation of the leading and next to leading infrared logs and the inclusion of the 
one loop corrections to the hard scattering amplitude (read one loop matching). In the fifth 
section we give our numerical results while the last section draws conclusions regarding what 
errors we can expect in the extraction process. 



2 Is Resummation Necessary? 

As mentioned above, the theoretical calculation of the lepton spectrum in inclusive decays 
breaks down near the endpoint. Both the non-perturbative as well as perturbative corrections 
become large in this region. Here we investigate the need to perform resummations in 
either or both of these expansions. The one loop decay spectrum including the leading 
non-perturbative corrections is given by |TD 



^1 - 



x^(3 - 2x)(l - + 2(3 - x)x^Eb - 1x^(9 + 2x)Ki, - \x^(lh + 2x)G'b 

Svr 3 3. 



+ 



E, - -K, + -G, 



5{l-x) + -Ki,5'{l-x). 



(1) 



Where 



31 5 2E 

lix) = log^(l — x) ^ log(l — x) + TT^ H — and x = — -. 

6 4 nih ' 



^0 — I Vub 



2 Glml 
' 967r3 ' 



(2) 
(3) 



2 



Eb, Gb and Kb are hadronic matrix elements of order and are given by 



Eb = Gb + Kb, 

Kb = {B{v)\K^^\B{v)), 

Gb = {B{v)\Kg'^^^f^K\B{v)), (4) 

Ami 

is the velocity dependent bottom quark field as defined in heavy quark effective field 
theory. From the above expressions we see that the breakdown of the expansions, in 
and e = ^^225. manifest themselves in the large logs and the derivative of delta functions, 
respectively. 

2.1 The non-perturbative expansion 

As one would expect for heavy meson decay, the leading order term in e reproduces the 
parton model result. All corrections due to the fact that the b quark is in a bound state are 
down by |@] . However, near the end point of the electron spectrum we begin to probe the 
non-perturbative physics. The general form of the expansion in e = to leading order in 
as, is given as follows 

^f- = e{l-x)(e' + e' + ---)+6{l-x)(0e + e' + e' + ---) (5) 
+ ■■■ + 6^"\l-x)(e^+' + e^+' + ■■■)+■■■ (6) 

The end point singularities are there because the true end point is determined by the 
mesonic mass and not the partonic mass, as enforced by the theta function in the leading 
order term. The difference between these end points will be on the order of a few hundred 
MeV. To make sense of this expansion we must smear the decay amplitude with some smooth 
function of x. Normally, this would not pose a problem, however, given that the distance 
between the b c and b ^ u end points is approximately 330 MeV , we are forced to 
integrate over a weighting function which has support in a relatively small region. On the 
other hand, if the weighting function is too narrow, then the expansion in e will not be well 
behaved. 
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Thus we must find a smearing function that minimizes the errors due to corrections 
which does not overlap with the energy region where we expect many b c events. The 
question then becomes how many b ^ c transitions can we allow without introducing large 
errors due to our ignorance of the b —>■ c end point spectrum (the theory breaks down in the 
b ^ c end point region as well though there are important differences between this case and 
the b —>■ u transitions)? 

The issue of smearing was addressed by Falk et. al.]^ who used Gaussian smearing 
functions to gain quantitative insight into the need for smearing. They found that without 
any resummation, the smearing function should have a width which is greater than e, but 
that after resumming the leading singularities, we need smear only over a region of width e. 
Here we will smear by taking moments of the electron energy spectrum (we work with the 
moments of the spectrum because it greatly facilitates the resummation of the perturbative 
corrections). Thus, we must address the question of what range of values of N will lead to 
a sensible expansion which is also not overly contaminated by 6 ^ c transitions? This will 
obviously depend on the ratio of Vub to He- To get a handle on the numerics, let us for the 
moment assume that we wish that the number oib ^ u transitions be at least equal to the 
number of 6 — > c transitions in our sample. In figure 1., we plot -^{N), which is the ratio 

he 

of mixing angles for which the A^th moments of the leading order rates for b ^ c and b ^ u 
transitions will be equal, and is given by 
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r™ x^ ^^'" 7, [6 - 3x„ + (x^ - 6)a; + IxV dx 
Jo (1 — x)'^ 
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Xm is for the B D transitions and takes the value x^ ~ .9. Given the bounds 
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.002 < < -024, (8) 

I Vbc r 

we see that an understanding of the spectrum for moments around ^ 20 is necessitated 
if we wish to keep the b ^ c contamination under control (Of course we do not suggest 
that these moments can be measured given the finite resolution of the experiment. We will 
discuss this situation later in the paper). 
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Figure 1: The ratio of for which the Nth moments of the leading order spectra are equaL 



We now consider the issue of determining the maximum value of for which the ex- 
pansion makes sense. Let us first consider the expansion in e. The moments of the leading 
singularities of eq.(^) will behave as 



TV 



Cn 



(9) 
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As a possible criterion on the size of A^, we may impose that there be no growth with n. 



That is 



Ar!e"+i (/ne + ^(A^-n + l)) 



< 0. 



(10) 



(A^-n)! 

e represents the value of some matrix element in the heavy quark effective theory. It is 
assumed that the value of e should be on the order of a few hundred MeV/nii,, but in theory 
it could vary by a factor of order one from term to term. To get a handle on the sizes of e, 
we may consider the leading e, which is given by 
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Quark model calculations suggest that e\ is on the order of .010. Thus, naively, it seems 
that to keep the expansion in e under control, we must keep A^ < 10. This estimate is perhaps 
too crude for our purposes given that we know nothing of the growth of the coefficients Cn 
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nor of the range of possible values of eA?, it does suggest that some sort of resummation may 
be necessary. 

Neubert0] pointed out that it is possible to resum the leading singularities, much as in 
the case of deep inelastic scattering, into a non-perturbative shape function 

f{k+) = {B\5{k+-tD+)\B). (12) 

This function gives the probability to find the b quark within the hadron with residual light 
cone momentum k+. Thus, this function is roughly determined by the kinetic energy of the b 
quark inside the meson. This structure function will be centered around zero and have some 
characteristic width 6. 6 will determine the maximum size of for which the expansion 
without resummation makes sense. To get a well behaved expansion we choose such that 

gives order one support to the structure function throughout its width. The value of 
6 is unknown at this time, and various authors have given different estimates for its value. 
We can assume that this width should be on the order of (m^ — mb)/2 which is around 300 
MeV. We shall choose, what we believe to be the conservative value of 500 MeV for 6. 
Since the structure function is the sum of derivatives of delta functions, we conclude that we 
should smear over the width of the function if we do not wish to incur large errors. Let us 
assume, for the sake of numerics, that should not fall below the value .1, within 500 MeV 
of the end point. Then we find that must be < 20. Thus, we expect the non-perturbative 
effects could be quite large for the range of that we consider here. Of course when A^ 
becomes very large, A^ > 100, it is necessary to go beyond leading twist since the soft gluon 
exchange in the t channel begins to dominate, not to mention the failure of the OPE due to 
its asymptotic nature p. 

We see that for our purposes we should include the non-perturbative structure function 
in our calculation. The fact that the knowledge of this non-perturbative function is needed 
to extract Vub should not bother us too much however, given that it is universal. That is 



to say we can remove it from our final result by taking the appropriate ratio ||TT|. Or it 
can be measured on the lattice, much in the same way that the moments of the proton 
structure functions are now being measured. Using the ACCMM model Blok and Mannel 
1^ concluded that a resummation of the non-perturbative corrections is unnecessary. If the 
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width of the structure function is smaller than the conservative number chosen here, then this 
could very well be true. This would be a welcomed simplification of the extraction process, 
since we would no longer need to rely on the extraction of non-perturbative parameters from 
other processes to measure the mixing angle Vub- 

2.2 The perturbative expansion 

Let us now address the issue of the perturbative corrections. The corrections in grow large 
near the electron energy end point, and, precisely at the end point, there are logarithmic 
infrared divergences. These divergences are due to the fact that near the end point gluon 
radiation is inhibited, and as a result, the usual cancelation of the infrared divergences 
between real and virtual gluon emission is nullified. Of course, the rate is not divergent, and 
we expect that a resummation procedure will have the effect of reducing the rate for the 
exclusive process. 

Near the end point large logs form a series of the form 

— = CuaLog'^ll — x] + Ci2aLog[l — x] + Ci^a 
ax 

+ C2ia'^Log'^[l - x] + C22a^Log^[l - x] H 

+ Cua^Log[l-xf + C^2a^Log^[l-x\^ 

+ . + . +■•• 

+ . + . +••• 

Which in terms of a moment expansion gives (for large A^)Q 

/ dxx^'^ = — (CuaLog'^N + CuaLogN + Ci^a 
Jo ax N ^ 

+ C2ia^Log^N + C22a^Log'^N H 

+ Csia^LogN^ + C^2a^Log^N +■■■). 

(13) 

^This form holds for b —>■ sj, for the scmi-leptonic decay we wih consider the moments of the derivative 
of the rate. 
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Figure 2: The difference between the moments given by the leading correction and the 
moments of the rate with only the double logs resummed. 

Given this expansion, we may ask what errors we expect to incur by truncating the 
expansion at order a^? For near 20, we see that 

—Log^N ~ .6, (14) 

TT 

SO we might expect that truncating at leading order would not be such a good idea. We 
must also note that in (|l|) the sub-leading log actually dominates the leading double log due 
to the large coefficient ^. The resummation of the double logs is simple and leads the the 
exponentiation of the double logs. Figure 2 shows the difference between the one loop result 
and the result with only the double logs resummed. We see that the difference is very small, 
on the order of five percent. Thus, one might come to the conclusion that no resummation 
of the perturbative series is necessary. However, given that the coefficients of the single 
logs as well as the vr^, which are just as large as the double logs for the range of we are 
considering here, are unknown at higher orders, we can only determine the errors induced 
by a truncation of the series after we have performed the resummation. 

Resumming the leading double logs in itself does not increase the range in over which 
perturbation theory is valid. Even after this resummation is performed the criteria for a 
convergent expansion is still ^Log'^N < 1 unless we know that the subleading logs expo- 
nentiate as well. However, one can show on very general grounds |]12| that all the end point 
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logs exponentiate as a consequence of the fact that these logs are really just UV logs in the 
effective field theory p4|, p]7| . Thus it is always possible to write down a differential equation 
for the rate based on its factorization scale independence. As such, the general form of the 
decay amplitude will be given by 



2n 

I fi .\ , ^ 

Log 

n=l m=l 



dx 



C(a,) + E E G^^Log^N. (15) 



Once we have this information, the question of the region of convergence becomes, do the 
lower order terms in question contribute numbers of order 1 in the exponent? We may 
continue to increase until we find that the subleading terms in the exponent contribute 
on the order 1. Thus in general resumming the leading logs does indeed allow us to take N 
into the range where ^Log'^N ~ 1. Here we will go further, as was done in [^, and sum 
the next to leading logs as well, allowing ^LogN < 1. This will allow us to determine the 
convergence of the expansion. Furthermore, we extend the analysis of P] to include the one 
loop matching corrections thus completing the calculation at order 

We wish to note that Blok and Mannel analyzed the effects of the large logs to the 
end point spectrum and concluded that no resummations were necessary. These authors 
propose to take the lower bound on the moment integral to be the charmed quark endpoint 
Xc- Doing this allows one to stay away from larger values of (the authors choose A^ < 10). 
Cutting off the integral introduces errors that have the doubly logarithmic Xc dependence 
/n^(l — Xc). To reduce these errors, it is necessary to go to higher values of A^. These authors 
claim that for A^ < 10, the errors induced by cutting off the integral are small, on the order 
of a few percent. However, we believe that these authors have underestimated their errors 
because they normalized their errors by the total width and not the moments themselves. 
Furthermore, and perhaps most importantly, the authors did not consider the possibility 
that the sub-leading logs could dominate the leading logs in the resummation, which as we 
shall see, is indeed the case. 

Finally, it should be pointed out that aside from being bounded by the size of the logs, 
A^ is bounded on purely logical grounds. The whole perturbative QCD framework loses 
meaning when the time scale for gluon emission becomes on the order of the hadronization 
time scale. This restriction bounds the minimum virtuality of the gluon, which we expect 
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to be on the order of ^(we will show this to be true when we perform the resummation) . 
Thus, performing resummations can only take one so far no matter how powerful one is. 
However, for top quark decays it is possible to get extremely close to the end point due to 
the large top quark mass. In this case it is clear that the resummation of the next to leading 
logs will become essential. Thus, the extraction of Vtd from inclusive top quark decays will 
have much smaller theoretical errors than in the b decay case. We shall discuss the issue of 
the breakdown of perturbation theory in greater detail after we perform the resummation. 



3 Factorization 

The large logs appearing in the perturbative expansion arise from the fact that at the edge of 
phase space gluon emission is suppressed. The problem of summing these large corrections 
has been treated previously for various applications, such as deep inelastic scattering and 
Drell-Yan processes|T^, 13], just to mention a few. The case of inclusive heavy quark decay 



has been treated previously in [Q]. An important ingredient of the resummation procedure 
is the proof of factorization. As applied to the present processes, this procedure separates 
the particular differential rate under consideration into sub-processes with disparate scales. 

In the case of inclusive semi-leptonic heavy quark decays, the relevant scales are and 
mb{l — x), with X = ^ in the rest frame of the b-quark. To understand how to best factorize 
the differential rate in the limit a; ^ 1, we need to know the momentum configurations which 
give leading contributions in that limit. With this in mind, let us consider the inclusive 
decay of the b-quark into an electron and neutrino of momenta Pe and Pi, respectively, and 
a hadronic jet of momenta ph- First we note that the kinematic analysis is simplified with 
the following choice of variables in the rest frame of the b quark pj 



x = — ?/o = y = 2 — ■ (16) 

nih 

The kinematic ranges for these variables are 

0<x<l; 0<y<x] {y /x + x) < yo < [y + I). (17) 
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Furthermore, define the variable 



1] = ( —] where x < rj < 1. (18) 

V2-Z/0/ 

This variable plays an analogous role to the Bjorken scaling variable in deep inelastic scat- 
tering phenomena. The invariant mass of the final state hadronic jet, and its energy are 
given by 

pl = mlil-r^)i2-yo), pl = '^{2-y,). (19) 

We should note that in determining the boundary values of the various variables we refer 
to the b-quark mass and not to that of the meson. This is justified within the perturbative 
framework we are working in at the moment. However, once we include the effects of the 
non-perturbative structure function, the phase space limits will take on their physical values. 

Let us now investigate the dominant momentum configurations near x ^ 1. First, we 
observe that the invariant mass of the hadronic jet + neutrino system is given by {p^ — Pe~ 
p^Y = ^bi^ ~ ^) which vanishes at the end point. The phase space configuration where the 
neutrino is soft is suppressed and hence, when the value of x approaches one, the electron 
and the hadronic jet-neutrino system move back to back in the rest frame of the b-quark. 
Furthermore, the invariant mass of the hadronic jet vanishes independently of the neutrino 
energy. This is readily verified using the phase space boundaries. The energy of the jet is 
large except near the point x ^ y ^ 1. In this region of the Dalitz plot factorization breaks 
down, and the techniques used here fail. However, this problematic region is irrelevant as 
a consequence of the fact that the rate to produce soft massless fermions are suppressed at 
the tree level. Thus, the following picture emerges at a; ~ 1. The b-quark decays into an 
electron moving back to back with the neutrino and a light-like hadronic jet. We choose 
the electron to be moving in the -|- (light cone) direction, and the jet moves in the — (light 
cone) direction in the rest frame of the b-quark. The constituents of the jet may interact 
via soft gluon radiation with each other and with the b-quark, but hard gluon exchange is 
disallowed. 



This simple picture is related by the Coleman-Norton theorem 0, |1^] to the type of 



Feynman diagrams that are infrared sensitive. According to this theorem, if we construct a 
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Figure 3: Reduced diagram for B decays. 



"reduced" diagram by contracting all off shell lines to a point, then at the infrared singular 
point, such a diagram describes a physically realizable process. Thus at x ~ 1 the type of 
diagrams that give large logs are precisely those described above and shown in figure ^. 

In the figure, S denotes a soft blob which interacts with the jet and the b-quark via soft 
lines. J denotes the hadronic jet and Hthe hard scattering amplitude. The typical momenta 
flowing through the hard sub-process are 0{mi,). H does not contain any large end point 
logs and has a well defined perturbative expansion in as{mf,). All the lines which constitute 
i/are off-shell and have been shrunken to a point. The soft function S contains typical soft 
momentum k , with k'^ k~ ^ k± ^ 0{mi,{l —x)). Thus, by "soft" we mean soft compared 
to m^, but still larger than Aqcd- The jet subprocess has typical momenta p such that 

^ with p'^ ■,p]_ ~ 0(mb(l — x)) and p^ ~ 0{mh). In order to delineate between 

momentum regimes, a factorization scale fi is introduced. The fact that the process is fi 
independent will be utilized to sum the large end point logs which are contained in the soft 
and jet functions. The reduced digram for the inclusive radiative decay b —>■ X^'j is exactly 
the same as above if we ignore the strange quark mass. 

An important consequence of the factorization is the fact that the soft function, S, is 
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universal. That is, it is independent of the final states as long as factorization holds. Thus 
the soft function in the semi-leptonic decay will be the same soft function as in the radiative 
decay. This universality will allow us to remove our ignorance of any nonperturbative physics 
due to bound state dynamics by taking the appropriate ratio. Thus, throughout this paper 
we will treat both the semi-leptonic as well as the radiative decays in turn. 

We conclude this discussion with a few comments. First, we should point out the differ- 
ences between factorization in the process considered here and in deep inelastic scattering 
for large values of the Bjorken scaling variable A crucial difference arises from the fact 
that the initial quark is massive, and hence, the semi-inclusive decays of the a heavy quark 
is infrared finite to all orders in perturbation theory because there are no collinear diver- 
gences arising from initial state radiation. This fact has the important consequence that the 
differential decay rate will be independent of /i. Whereas /i independence in deep inelastic 
scattering is only achieved after an appropriate subtraction is made with another process, 
such as Drell-Yan, which has the same collinear divergence structure as the deep inelastic 
scattering process. Next we note that, in general, the separation of diagrams into soft and jet 
subprocesses is not unique, and some prescription must be adopted. For a discussion of this 
issue see |12| . In our case, we will determine the proper separation from the requirement 
of the fi independence of the decay rate from the condition that the hard scattering ampli- 
tude does not contain any large end point logs, and that the purely collinear divergences in 
the jet must satisfy an Altarelli-Parisi like evolution equation. We will return to this point 
in the next section. The factorization can be made more manifest by going to the light like 
axial gauge with the gauge fixing vector pointing in the jet direction. In this gauge, the soft 
lines decouple from the jet on a diagram by diagram basis. 

In terms of the variables introduced earlier, the triply differential factorized decay am- 
plitude may be written as 



1 d^T r^inax 

Qm,{x-y){yQ-x) _ dk^f{k^ , fi')Jip^{pt - k+), fi')H{m,, fi'), (20) 



r° dxdydyo 



967r3 

This form will hold up to errors on the order of 0(1 — x 



•V^,\'ml (21) 
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We have chosen the electron to be travehng in the + direction with momenta k = 
mb{x, 0, 0_l), and 

f^max = ^bil - V), k^in = -{^3 - ^fe). (22) 

Here f{k~^) is the probabihty for the b quark to have hght cone residual momentum k^, and 
thus contains not only the information in the soft function S but also the non-perturbative 
information regarding the nature of the bound state. If we ignore perturbative "soft" gluon 
radiation, then this function coincides with f{k+) defined in the previous section. Notice 
that k"^^^ is negative. This is important non-perturbatively and represents the leakedge past 
the partonic end point due to the soft gluon getting energy and momentum from the light 
degrees of freedom inside the B meson. Loosely speaking it is due to the Fermi motion of 
the b- quark inside the meson. 

Less formally, we may write the derivative of the decay amplitude as 



-ld_ fdr\ /•2 
Fq dx \dx I 



j'dyo6{2-yoy{yo-l)G{x), (23) 



G{x)= dzf{z,mb/fi)J{mt{2-yo)iz-x),fi')H{mb{2-yo)/fi). (24) 

Jx 

In this equation we have changed variables from k^ to the residual light cone momentum 
fraction z = {1 — and absorbed a factor of ml into the jet factor. 

By taking the moments of this expression with respect to x we see that we are able to 
treat the hard, soft and jet functions separately. We are led to the following form for the 
moments of the semi-leptonic rate 

1 "'0 dx \dx J 

2 

dyo6{2 - yo){yo - l)fNJN{ml{2 - y^), ^i^)H{mi,{2 - yo),fi), (25) 

f Mb /nib 

In = / z'^dzfiz), (26) 

Jo 

Jjv(2-i/o)= f' X''j{ml{2~yo)il-X,fi'))dX. (27) 
Jo 
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In writing the last few equations we have dropped all terms of order 0(1 — x) or equivalently 
taken the large limit. The left hand side of Eq(p5D defines the moments of the semi-leptonic 
decay electron distribution, M§. 

The moments of the soft function may be decomposed into a product of moments 
of a perturbatively calculable ctat and the non-perturbative structure function Sn, which 
corresponds to the moments of / discussed in the previous section. We may write 

m = / -S{yM-), (28) 

and thus, taking moments, 

In = ctnSn- (29) 

An analogous situation exists for the decay B Xg'y. We define x = in the rest frame 
of the b-quark, and take the photon to be moving in the + direction, and, as in the semi- 
leptonic decay, at x ~ 1 the hadronic jet is moving in the — direction. Furthermore, the 
invariant mass of the hadronic jet and its energy are 

pl = ml{l~x), p° = m,(l-x/2). (30) 

Thus the s-quark is very energetic and since the invariant mass of the jet vanishes as a; ~ 1, 

2 

the s-quark decays into quanta which are coUinear once we ignore effects on the order of 

h 

Clearly the factorization picture discussed earlier for the semi-leptonic decay holds here as 
well and the reduced diagram is the same as in fig(3). 

As before we may take the moments of the differential rate 

Ml^— rfxx^-i— = SmOnJn- (31) 

1 ^ JO ax 



where 31 



= I V,,V,: p C',{m,). (32) 



and 

fi 

.N-l T/_2/i „ \ ,.2 



Jn= / dyy''-'J{mi{l~y),ii'). (33) 
Jo 

(Tat and Sn are the same functions defined in (^) and C-j is the Wilson coefficient of O7 as 
defined in For the radiative decay, the ln{\ — x)/(l — x)+ distribution in the amplitude 
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will correspond to ln^{N) in the moment. Whereas, in the semi-leptonic decay, taking the 
derivative of the amplitude will generate plus distributions which will then generate LogN 
and Log^N after taking the moments. Thus, we have reduced the problem of the resumma- 
tion of the large logs in the amplitude to resumming the logs in and ctat separately. This 
greatly simplifies the calculation as will be seen below. 



4 Resummation 

The resummation of the infrared logs is analogous to summing ultra violet logs. One takes 
advantage of the /i independence of the amplitude. In the case of infra-red logs, /i is the 
factorization scale, or equivalently the renormalization scale within the appropriate effective 
field theory, which for this case would the field theory of Wilson lines [P^. 

We first outline the derivation of a representation of the soft function near x = 1 
following the techniques developed in reference [^. We will work in the eikonal approxi- 
mation where soft momenta are ignored wherever possible. At the one loop level, the real 
gluon emission contribution factorizes and the quantity multiplying the tree level rate is 



The 5-function enforces the phase space constraint. Similarly the one loop virtual gluon 
contribution is given by 

F-i^ = -9'C,S{1 - X) [ - -4t^| . (35) 



. (27r)32A;o [p-kq-k {p ■ kf ^ 

Where p and q are the b quark and light-quark momenta respectively. In the Abelian 
theory exponentiation follows simply as a consequence of the factorization in the eikonal 
approximation. For each gluon emission one gets a factor of F™J1^ which is unitarized by the 
virtual contribution. After appropriate symmetrization the exponentiation follows. Next 
we use the result that even in a non-abelian theory, for the semi-inclusive process under 



consideration, exponentiation of the one loop result takes place |23, 24|. By considering the 
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jyth ]2ioment of the soft part, we obtain 



CTtv 



2p ■ q 



(27r)32A;o 



mt 



2ko 
M 



(36) 



p ■ k q ■ k (p ■ ky 

It should be noted that the ultraviolet cutoff is determined by the factorization scale 
fi. This cutoff is necessary despite the fact that the process under consideration is infrared 
finite. All momentum above this scale get shuffled into the hard scattering amplitude H. 
The need for a cut off stems from the fact that we have used the eikonal approximation. 
This approximation is equivalent to a Wilson line formulation of the problem, and thus, as in 
heavy quark effective field theory, generates a new velocity dependent anomalous dimension 

By an appropriate change of variables may be written as 



(Tat 



exp 



rfScA{kl)-C,^{mW))}. 

Jm?v^ k, n TT J 



y 

a 
n 



(37) 



In arriving at the above, we have made the replacement a.^ —>■ asikjj. This change has the 
effect of resumming the next to leading logs coming from collinear emission of light fermion 
pairsp6|. However, it does not sum all the soft sub-leading logs. 

Explicit calculations carried out at the two loop level |]T^ ^ indicate that the rest of 
the sub-leading terms in the above may be included |]13 

Cpasikl) 



TT 



Aiasiki)), 



(38) 



with 



A{as 



TT 



—Cf + — -Cpk 
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18 



TT 

~6 



10 



9 



TrNj. 



(39) 



This resums all the leading and next to leading logs of N in the soft function. 
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Thus, we obtain 



aNimb/fi) = exp{- 



y 



ii-ii-yr-') 



'^my kl 



■A{as{kl)) + B{as{my))}, 



(40) 



with B{as) = —as/iT. This integral is not well defined due to the existence of the Landau 
pole, and a prescription is needed to define the integral. Choosing a prescription leaves an 
ambiguity on the order of the power corrections 0|2y][^. If we use the large N identity 

1 



„N-1 



X 



1 



X 



N 



N 



N 
iVo' 



iVo 



(41) 



which is accurate to within 2% at = 10, to rewrite 

dk 



exp< 



rUb/N kj 



2A{a,{k^))ln'^ - B (asikl)) 



(42) 



then we have fixed a prescription which is unambiguous to the accuracy we are concerned 
with in this paper. From this result, we find that aN{mb/fJ^) satisfies the RG equation 



J^dfi dgl \ , 



2A (asii/))\n^ + B(asifJ.' 

^ ' VTLb ^ 



(43) 



We note that this is in agreement with reference]^] where Wilson line techniques were utilized. 

We may now use the /i dependence of the soft function, together with the fact that the 
total amplitude is /i independent, to determine the renormalization group equation satisfied 
by the jet and hard functions. We have seen in section (3) that the A^*^ of the derivative of 
the moments of the semi-leptonic decay has the factorized form 



(JN{mb^2 - vqI ^i)jfj{mb/ fi)H'\ml{2 - yo)/^^). 



(44) 



Whereas, for the radiative decay the moments of the decay spectrum is given by 

aN{mb/fi)Jjf{mb/fi)H'^{ml/fi'^). (45) 
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We have now labeled the jet and hard functions according to their processes since these 
function are not universal. We will first consider the RG equation satisfied by Jjj and . 
The equations satisfied by J^^ and i?*' can then be determined by simply by making the 
appropriate replacements. We may derive the RG equations satisfied by these functions by 
using the following facts 

• lJ'-^{(^N{mb/ li)JN{mb/ fi)H{mb/ fi)) = 

• The RG equation satisfied by aj\f{mb/fi) is given by Eq. (|43| ) 

• The hard scattering amplitude by definition has no dependence 

• The jet functional form which is Jn 
Leading to the following RG equations 



N 



( d d \ li^N 

yu— + (3— + f{as) Jn {nib/fi) = 2A(a,(/i^)) In — ^JNinib/n), (46) 

^ ~ ^^"'^ " B{a,{i2'))^ H{mb/fi) = -2A (^a,{fi') In H{mb/fi). (47) 

f{as) is an arbitrary function which can only be determined from additional input. We fix 
f{as) by requiring that the purely collinear divergences of the jet factor be determined by 



an Altarelli-Parisi type equation as discussed in [|T3|, [T^. We note that for these purposes 
the jet factor is a cut light quark propagator in the axial gauge. 

By requiring that we correctly reproduce the pure collinear divergences at the one-loop 
level it is found that 

/(«,) = 27(«,). (48) 
Where jic^s) is the axial gauge anomalous dimension []T^, ^ 



3 rv 

7(«.) = -t-Cf + .... (49) 

4 vr 

The solution of the jet RG equation may be written (to the desired accuracy) in the form 



Jjf{mb/fi) = exp 



iv K 



(50) 



19 



For future purposes, we rewrite this in the form 

o I o 1 r 



Jl{mb/ii) = expi / — 
\Jo y - 



y 



dy" 



A{a,{kl)) - ^{a,{mly)) 



■A{a,{kl)) - j{a,{mly)) 



(51) 



We may now write the exphcit expressions for the resummed jet and soft factors. For 
the radiative decay B Xgj we rewrite the various representations obtained earher leaving 

'■I dz 



Mrrib/fi) = exp / -^(1 - z''-')N{z), 
Jo 1 — z 



(52) 



JN{mb/ii) = exp [ -^(1 - z^ ^)I{z), 
Jo i — z 



(53) 



N{z) = 



r 

Jml(l-z)^ 



dk\ 



kl 



A{as{kl)) 



B{as{ml{l - zf)) 



j-H/mb rn~ 
Jl J mly' 



A{as{kl)) + B{as{mly^)) 



(54) 



I{z) 



Jmi 



de, 



mlil-z) kj 



+ 



ri^ dy 
Jl y 



A{a,{kl))]-^{as{mUl-z))) 
dk^ 



k' 



■-A{as{ki))-'y{a,{miy)) 



(55) 



Combining these two factors we see that for the N dependent piece in the exponent, 
the /i^ dependence exactly cancels. There are, however, pieces which are independent of 
N which are /i dependent and these will combine with similar terms in the hard scattering 
amplitude to give a /i-independent answer which must be true by construction. Combining 

all the factors we find 



aN{mb/ii)JN{mb/ii) = exp - / (1 - z^~^)K{z), 

Jo 1 — z 



(56) 
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K{z)^ -^A{as{kl))-^{as{ml{l-z)))-B{as{ml{l-z)%. (57) 

The N^^ moment of the decay rate in then given by 

Ml^SMCJNJjjH\as{ml)). (58) 

The value of the the one loop hard scattering amplitude H'^ is given in the Appendix. 

For the case of the semi-leptonic decay the expression for the soft factor is the same as 
above. However, for the jet, we must rescale J^{mb/ ji) Jjv {m},/ 11^/2 — yo) ■ Thus, we get 

(m,/^i^/l^) = expl^f^ ^)^-i] L(y,a;.)|, (59) 

L{y, x.)^ ^A{a,{kl)) - 7(«.(m^y(l - x^))). (60) 

In writing the above, we have used the fact that = x+x^, and a; ~ 1 to replace the variable 
I/O by Xu, the neutrino energy fraction. After some algebra the above may be combined with 
the expression for the perturbative soft function, such that for the product we may write 

<7^(— )J^(— yr^)=exp|- f'^(l-z''-')Q{z) + P{N,xA, (61) 
[1 n y JQ i — z J 

Q{z) = / -TT^Mkl)) - 7(«.(m^(l - z))) - B{a,{ml{l - zf)\ (62) 



^2(1-^)2 ki 



7 f«s (-^j) -^iMmly)) 



Ji/N y Jmlv{i-x^) ki Ji-x. y 

(63) 

We note that in deriving this form we have kept only the N dependent pieces in the 
exponent. Our analysis shows that certain N independent terms, like those proportional 
to ln(l — Xjy), can also be resummed using the above mentioned procedure. However, we 
have taken as{ml) ln(l — Xj,) to be small in the relevant Xi, range and hence relegated all of 
these logs to the hard scattering amplitude. Thus, we may write for the N^'^ moment of the 
semi-leptonic decay rate, up to corrections 0{1/N), as 



pi 

^ ^ Sn dx„6x„(l - Xt,)aNJN(ml(l - x^,), f/)H(mb(l - xJ, /x). (64) 
Jo 
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The complete expression for the one loop hard scattering amplitude both the radiative and 
semi-leptonic processes at x ~ 1 are given in the Appendix. 

We conclude this section by giving some simplified expressions for the product (TnJn 
which will be useful for numerical analysis. We begin by noticing that, as long as as{ml)lnN < 
1, the resummation formulae given above have a convergent power series expansion in 
as{ml)lnN to the next to leading log accuracy. Thus, we compute these expressions to 
this accuracy and delegate all the non perturbative effects phenomenologically to the struc- 



ture function Sn- For a similar approach for the case of e~^e~ annihilation see fl^. To 
evaluate the integrals in the exponent, we may perform the z integration using the large N 
identity (^T]). The k± integration is simplified by using the RG equation for the running 
coupling to change variables to as, i.e. 

^-^ = -T^-^(^-l--^ + Oia^s)) (65) 



where 



^ llCA-2Nf ^ nCl-bCANf-SCpNf 



Next we use the expansion, correct to next to leading log accuracy, 

c^s{ml/N) = "-^7'^ f 1 - ^ ""'^fl Inil - PoasimDlnN)] , (67) 
1 - l3oas{mf)lnN \ Pol- Poas{mi)lnN J 

to obtain 

Har,JN) = InN {gj{x)) + gl{x), (68) 

where, 

X = f3oas{ml)lnN. (69) 

In the above, the functions gi and g2 have the following form for the two processes discussed 
in this paper 

For the radiative decay 

9l = - 2x)/n(l - 2x) - 2(1 - x)/n(l - x)), (70) 
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and 

92 = -^(-Mi-2x) + 2Mi-x))-4^(Mi-2x)-2Mi-x) 

+ \ln\l - 2x) - ln\l - x)) + ^-^H^ - x) + " 2x) 

2 ^ vr/^o 27r/^o 

^(1) 

- -—lnN^{ln{l-2x)-ln(l-x)). (71) 
For the semi-leptonic decay 

gf = 9l (72) 

and 

92 = 9l + - Xu)Hl - x)- (73) 

We have kept only the N dependent terms in these g factors which exponentiate. There are 
also Nq dependent constant terms which we will shuffle into the hard scattering amplitude. 
These terms are given by 

h'^ = + A^'^^lnNoln{l-x,). (74) 

TT 

Furthermore eq. (|6^ ) for the semi-leptonic case becomes 

- — / x""-^ — — = Sn / dx,Q{l - x,)x,{H{x,) + h'')Exp{gi + g^). (75) 
i "'0 dx dx Jo 

In writing the above, we have used the notation 

A(a,) = (^)A« + (^)2A(2), (76) 

71 71 

and 

Bia,) = (^)i?«, 7(«s) = (-)7(^)- (77) 

The values for A^^\A^^\ 5^^), 7(1) have been given previously. H[xy) is the one loop correc- 
tion to the hard scattering amplitude given in the Appendix. 

It is interesting to note that if we expand the expressions for gf and gf in (0), we 
see that G24, as defined in (|l^), vanishes. Thus, the two loop results does not trivially 
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exponentiate as one might have naively thought. Such behavior is a universal property of 
the asymptotic limit of distribution functions and is a consequence of the fact that only 
"maximally non-abelian" graphs contribute to the exponent beyond one loopp^. Knowing 



this greatly reduces the number of graphs that need to be calculated in a general resummation 
procedure. 

From expression ( |75D we may determine the range of for which our calculation is valid. 
The integration over y contains a branch cut at = ^^^^ , signaling the breakdown of the 
the perturbative formalism. This breakdown is coming from the fact that the time scale for 
gluon emission is becoming too long. An inspection of the resummation formulae for these 
quantities suggests that in the region z ~ 1 such that fc^ < A^ , non-perturbative effects 
become important. Thus, we conclude that we may only trust our results in the range 

N < (78) 

^^QCD 

5 Analysis and Results 

With the resummation now in hand, let us consider the relative sizes of all the contributions. 
In figure § we show the difference between the one loop result and the resummed rate given 
by eq.([75|) normalized to the moments of the one loop result, 

- — / x""-' — — = 1 - TT^ + - - -InN + In^N). 79 
To Jo dx dx Stt 4 6 

In our calculation we take Ag^T^ ^ 200 MeV. We see that resumming the next to leading logs 

has a 20% — 50% effect in the range of A^ we are considering. Furthermore, for completeness 

we have included the effect of resumming the vr^. The result of this resummation is given 

by the dashed line. We see that the effect of resumming the vr^ is smallQ. As a check of the 

numerics we compared the resummed expression to the one loop result ( [79D and found that 

for small A^ the two coincide to within less than a percent. The fact that the resummation 

of the next to leading logs is more important that the leading logs, is rather disheartening. 

^Note that in resumming the tt^ we only resum part of the tt^ in the expression (2), since part of the tt^ 
contribution comes from integration over the neutrino energy 
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Figure 4: The difference between the moments of the one loop result and the resummed 
result with (dashed) and without (undashed) the resummation of the tt^, normalized to the 
one loop result. N varies from 10 to 30. 

It leads one to believe that perhaps the next to next to leading logs will be even more 
important. However, the fact that the effect of subleading logs is larger than the leading is 
already hinted at one loop, given that the ratio of the coefficients in front of these logs is 
31/6. It could be hoped that the ratio of the coefficients of the next to leading and next to 
next to leading logs is not so large and the terms left over in our resummation will be on the 
order of 10%. 

6 Discussion 

Before we conclude with a discussion of the future prospects of the extraction of Vub we wish 
to point out that there is one tacit assumption which has been made up to this point in 
our investigation. That is, we have assumed that local duality will hold when we are a few 
hundred MeV from the end point. The whole formalism of using the OPE in calculating 
inclusive decay rates assumes that at certain parts of the Dalitz plot, the Minkowski space 
calculation will give the correct result. This should be a good approximation as long as we 
stay away from the resonance region. The question is, how far from the end point does this 
region begin? If it is found that single resonances dominate, even as far as a few hundred 
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MeV from the endpoint, then the extraction of Vub through inclusive decays is surely doomed. 
The quark model seems to indicate that this may be the case though other theoretical 
predictions say otherwise. We will have to wait to see the data before we can decide on the 
fate of the extraction methods discussed here. 

Next we wish to reiterate that completely eliminating the background from b ^ c transi- 
tion by going to very large > 30, is not feasible since there is no way to reliably calculate 
the soft gluon emission which takes place. This is because when one goes that far out on 
the tail, the time scale for gluon emission is too long compared to the QCD scale to have 
any hope of perturbation theory making any sense. Again, this statement is independent 
of how many soft logs one is willing to resum. Another way of saying this is that when x 
gets to close to one, there is no operator product expansion since the expansion parameters 
is • Thus, we are stuck with the fact that there will always be contamination from 

transitions to charmed final states. Calculating the end point of the charmed spectrum using 
the techniques discussed above fail as well since resonances will dominate. Thus, there does 
not seem to be any way to avoid having to use a model to determine the background in the 
extraction process. The best we can hope to do, using the results in this paper, is to go to 
a large enough value of A^ that we can reduce the model dependence as much as possible. 
Certainly, we can greatly reduce the model dependence from what it is in present extractions 
which rely solely on models. 

The last point that needs to mentioned is the fact that measuring large moments it not 
experimentally feasible, as varies much too rapidly. For instance, if we assume that the 
bin size is given by 5 = then the error at point x for the A^th moment will be 
Therefore, the error can accumulate quite rapidly. Thus, it will be necessary to take the 
Mellin transform of our result. Given that our result is only trustable for A^ < 25, one must 
be careful to calculate the contribution to the inverse transform from higher moments, if one 
hopes to impose the bounds on the errors discussed in this paper. Also, for smaller values 
of A^ one must be sure not to use the resummed formula as we have dropped terms that go 
like 1 — X. 

Given these caveats, we may now address the issue as to what accuracy we can deter- 
mine Vub using inclusive decays. Since the sub-leading logs dominate the leading logs, the 
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conservative conclusion would be that a model independent extraction of Vub is not possi- 
ble. However, let us proceed under the assumption the sub-sub-leading logs will be smaller 
than the sub-leading logs. In this case we may say that we have been able to reduce the 
errors from radiative corrections down to the order of 10%. However, the QCD perturbative 
expansion is notoriously asymptotic, and though we may hope that we have resummed the 
dominant pieces of the expansion, there could still be large constants (independent of A^) 
which could arise. 

Another source of errors will come from the fact that we need to eliminate the dependence 
of the decay rate on the moments of the non-perturbative structure function by taking the 
ratio of the semi-leptonic decay moments with the moments of the radiative decay. This will 
introduce the errors in the radiative decay into the semi-leptonic decay. One could calculate 
without any non-perturbative resummation, thus eliminating these errors (the results in this 
paper are easily modified to include this possibility), but then it is difficult to quantify the 
model dependent errors introduced in the truncation. Finally, there are the errors introduced 
due to the model dependence from the calculation of the background. This error will be 
reduced as we choose larger values of A^. This is the most difficult error to quantify, and we 
shall not discuss it here. 

The authors believe that, if the end point is not dominated by single particle resonances, 
and if we assume that the fact that the sub- leading logs dominate the leading logs is just 
an anomaly, then we may hope to eventually extract Vub at the 30% level using the results 
presented here. Moreover, resumming the next to leading logs is indeed necessary. However, 
the more conservative view would be that the endpoint calculation is just intractable at 
this time , since it could be that the sub-sub-leading logs will dominate. To be sure that 
this is not the case the sub-sub-leading logs would need to be resummed. This would entail 
calculating A to three loops, and B and 7 to two loops. Without this calculation, we can 
not determine with certainty the size of the errors. 
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7 Appendix 



In this Appendix we give explicit expressions for the hard scattering amphtude at the one 
loop level and to leading order in (1 — x). We first present the results of the computation 
of the QCD corrections to the doubly differential rate J^J^^ for the semi-leptonic b-quark 
decay. It is clear from sections 2,3 that this is the quantity whose moments factorize, and 
which is relevant for the resummation. We write, 

-f — = 6ro(2/o - 1)(2 - 2/o)(l - -^G{x,yo)). (80) 
axayo Sir 

where, Tq was defined earlier. The contributions of the real and the virtual gluon emission 
diagrams to G{x, yo) are given by 

G}Tn = In^il + - x) - 2/n(2 - yo)ln{l - x) - ln\2 - y^) + in(2 - Z/o) + ^ 

■'2 2d 

\ln\^) - ^/n(4) + 2/n(2 - yo)ln{^), (81) 
2 2 mi 



and 



Gj: = ^ln\2-y,)-2lni2-y,)ln{y,-l) + 2Re{Li2{- )) + 3^ ^/n(2 - yo) 

Jn{2-yQ) 5 tt^ 1, 2/A\ 5, ^A^, 
?/o — 1 2 6 2 mi 2 mf 

- 2ln{2-yo)ln{^). (82) 

In the above, A is the gluon mass used to regulate the infrared divergences at intermediate 
stages of the calculation. Combining these results and integrating over yQ gives the electron 
spectrum which agrees with but disagrees with 

From this we see that to the approximation we are working in, the hard scattering 



amplitude as defined in eq.^ is given by 

Hsi = l-^(^^ + 2ln\2-yo)-2ln{2~yo)ln{yo-l)-hn{2-yo) 

, ^.2.,.Mj^)).a (S3) 
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For the radiative decay, we may extract the hard scattering amphtude from 

2a, ,3 27r\ 
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